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Abstract While tested to a high level of accuracy in the Solar system, general 
relativity is under the spotlight of both theoreticians and observers on larger 
scales, mainly because of the need to introduce dark matter and dark energy in 
u: the cosmological model. This text reviews the main tests of general relativity 

focusing on the large scale structure and more particularly weak lensing. The 
!■ complementarity with other tests (including those on Solar system scales and 

the equivalence principle) is discussed. 
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1 Introduction 



Gravitational lensing is historically bound to the developments of general rela- 
| tivity (GR) and, more generally, of the theories of gravitation. Since the end of 

. the 18th century, it was thought that light can be deflected by a gravitational 

field, in particular with the works of Georg von Soldner that postulated that 
qq . light must behave as any other particle or of Robert Blair, John Mitchell and 

Pierre Simon de Laplace (see Ref. p] for an historical discussion). 
<3\ . The deflection of light by any massive body is a central prediction of GR. 

In particular, the observations of the deflection of light emitted by distant stars 
by the Sun during the Solar eclipse on the 29th May 1919 by the expeditions 
led by Eddington and Cottingham on Principe island and by Davidson and 
Crommclin in the Nordeste region of Brasil is always considered as an exper- 



imcntal confirmation of the predictions of GR. Indeed, if such an observation 
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was a test of GR, it is only because the mass of the Sun was supposed to 
be well-determined at the time. On the one hand, the light deflection angle 
predicted by GR is 

_ 4GM 
- ~b^~ 

where G is the Newton constant, M© the Solar mass, b the impact parameter 
and c the speed of light, while, on the other hand, the dynamics of massive 
bodies, such as planets, are in an extremely good approximation still given by 
Kepler third law, 

P a? 



2tt GM & ' 

where P is the period of the orbit and a its semi-major axis. Measuring A9 
and b on one side and P and a on the other gives two estimates of GM© , 

ions _ be 2 AO dyn _ 27ra 3 



GM' cns = , GM, 



4 > w p 

that needs, given the error bars, to be consistent. This illustrates that lensing 
alone does not allow to construct a test of the theory of gravity but that we 
need to check the consistency between various predictions such as lensing and 
dynamics of massive bodies. 

Today, gravity, i.e. the only long range force that cannot be screened, is de- 
scribed by GR in which it is the consequence of the geometry of the spacetime. 
General relativity is consistent with all precision experimental tests available 
but most of these tests are restricted to the Solar system or to our Galaxy, 
so that they are only local. By considering astrophysical systems, we can ex- 
tend the domain of validity of GR at large distance, low acceleration or low 
curvature. In particular, most atttempts to construct a quantum theory of 
gravity or to unify it with other interactions predict the existence of partners 
to the graviton, i.e. extra fields contributing to a long range force, and thus 
to gravity (e.g. this this the case in all extra-dimensional theories where some 
components of the extra-dimensional part of the metric behave as scalar fields 
from a 4-dimensional point of view; string theory also predicts the existence 
of a scalar field, the dilaton, in the graviton supermultiplet). It follows that 
deviation from GR [2] are expected but in many cases (such as scalar-tensor 
theories) , the theory can be dynamically attracted toward GR so that cosmol- 
ogy can set sharper constraints than those obtained locally. 

An other reason to test GR in these regimes is related to our current cos- 
mological model. This model, which is consistent with almost all astrophysical 
data requires the addition of both dark matter (a fluid with negligible pressure 
that does not interact with standard matter) and dark energy (a fluid with a 
negative pressure), which represent, respectively 23% and 72% of the matter 
content of the universe. The need for these two components arises from the 
study of the dynamics of clusters, galaxies, large scale structures and of the 
cosmic expansion under the assumption that GR holds on astrophysical scales. 
This conclusion has been challenged by invoking possible modifications of GR 
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either in a low acceleration regime to explain the galaxy rotation curves and 
cluster dynamics without introducing dark matter and at large distance to 
account for the late time acceleration of the cosmic expansion. 

In conclusion, the possibility to sharpen our understanding of the validity 
of GR from astrophysical data and the need to understand the properties of 
dark energy and dark matter, which are tied to the validity of GR, are our 
two driving motivations. The main difficulty is that, on astrophysical scales, 
most observations entangle the properties of gravity, matter, as well as other 
hypothesis such as the Copernican principle. 

The bottom line of the construction of these tests is simple. Once GR is 
assumed valid, it describes the dynamics of the cosmic expansion, the growth of 
large scale structures, the propagation of light etc... so that many observables 
are not independent. Such observables can be used to construct consistency 
relations. Any sign of a violation of these relations will indicate the need to 
extend our description of gravity, but will not indicate us how. For instance, 
in the oversimplified Solar system example described above, we want the two 
estimates of GMq to agree that is we must have 

be 2 A9 27ra 3 _ 
4 P 

This is a relation between observable quantities (b, A9, P, a) that has to be 
satisfied. Actually, it was first proposed in Ref. [3J to perform a similar test 
on cosmological scales using weak lensing and galaxy redshift suveys, followed 
by the analysis of Ref. [3] . 

The review is organized as follows. We start, in § [2j by recalling the main 
hypothesis on which GR is based as well as the standard constraints obtained 
in the Solar system. We also discuss the use of alternative theories and draw 
the conclusions of what was learnt in the Solar system for constructing tests 
on astrophysical scales. In § [31 we discuss briefly tests on galactic and cluster 
scales where the need for dark matter can be interpreted as the necessity to 
modify GR in a low acceleration regime. Larger scales are considered in § 0] 
which focuses on the large scale structure of the universe. 



2 Relativity and its Solar system tests 

2.1 General relativity (in brief) 

Let us recall that GR, Einstein's theory of gravity, relies on two independent 
hypothesis. 

First, the theory rests on the Einstein equivalence principle, which includes 
the universality of free fall, the local position and local Lorentz invariances in 
its weak form (as other metric theories) and is conjectured to satisfy it in its 
strong form. We refer to Ref. [5] for a detailed description of these principles 
and their implications. The weak equivalence principle can be mathematically 
implemented by assuming that all matter fields, including gauge bosons, are 



1 



minimally coupled to a single metric tensor g M „ . This metric defines the lengths 
and times measured by laboratory rods and clocks so that it can be called the 
physical metric. This implies that the action for any matter field, tp say, is of 
the form 

S m [i>,9nu]- (2) 

This so-called metric coupling ensures in particular the validity of the univer- 
sality of free-fall. 

Then, the action for the gravitational sector is given by the Einstein-Hilbert 
action 

c 3 



"Sgravity — ^q^q J ^ X ^J 9*R*i (3) 

where g* is a massless spin-2 field called the Einstein metric. The second 
hypothesis of GR states that both metrics coincide, i.e. 

* 

9y,v — 9fiv 



2.2 Testing GR 



It follows that one can aim at testing both the equivalence principle and the 
dynamical equations that derive from the Einstein-Hilbert action. 

The assumption of metric coupling is well tested in the Solar system. First it 
implies that all non-gravitational constants are spacetime independent, which 
have been tested to a very high accuracy in many physical systems and for 
various fundamental constants 6,7,8,9 , e.g. at the 10 -7 level for the fine 
structure constant on time scales ranging to 2-4 Gyrs. Second, the isotropy 
has been tested from the constraint on the possible quadrupolar shift of nu- 
clear energy levels [rUlfTTUT^] proving that matter couples to a unique metric 
tensor at the 10 -27 level. Third, the universality of free fall of test bodies in 
an external gravitational field at the 10 -13 level in the laboratory [13, 14 . The 
Lunar Laser ranging experiment |15j , which compares the relative acceleration 
of the Earth and Moon in the gravitational field of the Sun, also probe the 
strong equivalence principle at the 10 -4 level. Fourth, the Einstein effect (or 
gravitational redshift) states that two identical clocks located at two different 
positions in a static Newton potential U and compared by means of electro- 
magnetic signals shall exhibit a difference in clock rates of 1 + \Ui — E^/c 2 . 
This effect has been measured at the 2 x 10 -4 level [IB] . 

The parameterized post-Newtonian formalism (PPN) is a general formal- 
ism that introduces 10 phenomcnological parameters to describe any possi- 
ble deviation from GR at the first post-Newtonian order [5]. The formalism 
assumes that gravity is described by a metric and that it does not involve 
any characteristic scale. In its simplest form, it reduces to the two Edding- 
ton parameters entering the metric of the Schwartzschild metric in isotropic 
coordinates 

2Gm PPN ( 2Gm \ 2 ( PPN 2Gm\ 
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Indeed, general relativity predicts /3 = 7 PPN = 1. These two phenomeno- 
logical parameters are constrained (1) by the shift of the Mercury perihe- 
lion [17] which implies that |2 7 PPN - /3 PPN - 1| < 3 x 1CT 3 , (2) the Lunar laser 
ranging experiments [T5] which implies |4/3 PPN — 7 PPN — 3| = (4.4±4.5) x 10~ 4 
and (3) by the deflection of electromagnetic signals which are all controlled 
by 7 PPN . For instance the very long baseline interferometry [18) implies that 
|7 PPN — 1| = 4 x 10~ 4 while the measurement of the time delay variation to 
the Cassini spacecraft [19] sets 7 PPN - 1 = (2.1 ± 2.3) x 1(T 5 . 

The PPN formalism does not allow to test finite range effects that could 
be caused e.g. by a massive degree of freedom. In that case one expects a 
Yukawa-type deviation from the Newton potential, 

l + a(l-e- r / A )] , 

that can be probed by "fifth force" experimental searches. A characterizes 
the range of the Yukawa deviation while its strength a may also include 
a composition-dependence [5J. The constraints on (A, a) are summarized in 
Ref. [20] which typically shows that a < 10~ 2 on scales ranging from the 
millimeter to the Solar system size. 

GR is also tested with pulsars [2T1I22] and in the strong field regime |23) . 
For more details we refer to Refs. 5,24,25]. Needless to say that any extension 
of GR has to pass these constraints. However, these deviations can be larger in 
the past, as we shall see, which makes cosmology an interesting field to extend 
these constraints. 




2.3 Alternative theories of gravity 

The ways of modifying GR are so various and in large number that we cannot 
review them here. We refer to Refs. [5,26 for some examples. 

Let us however introduce the simplest modification of GR, that is scalar- 
tensor theories of gravity in which gravity is mediated not only by a massless 
spin-2 graviton but also by a spin-0 scalar field that couples universally to 
matter fields (this ensures the universality of free fall). Their action can be 
written as, in the so-called Einstein frame, 

S = -^~c~ J d 4 x^T4R* - 2gF'd„<p*d v <p* - W(<p*)] 

+ Smattcr[^ 2 (^*)s;L;# (4) 

where G* is the bare gravitational constant. The physical metric, to which 
matter is universally coupled, <? M „ = A 2 {ip if )g*^ v is the product of the coupling 
function A, which characterizes the strength of the scalar interaction, and 
the Einstein frame metric This theory involves a new degree of freedom 
coupled to matter. 
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It can be used to illustrate the effect of modification of GR on lensing 
(see Ref. [26] for more details). Consider the action for electromagnetism in d 
dimensions 



^Maxwell = T / V~9*9* V 9^" FupF va d d X 



transforms to 



^Maxwcii = \ J V=g~g> u 'g pa A d -*F tip F va d d x 

under the confomal transformation g^ v — A 2 (ip*)g^ L ,. In the relevant case of 
a d — 4 dimensional spacetime, the Maxwell action is conformally invariant. 
Therefore light is only coupled to the spin-2 field g* so that light deflection 
by a point mass M must be the same as in GR, i.e. 

An AG* MA 2 
A9 = 



be 2 ' 

where A 2 M is the deflecting mass in the Einstein frame. It thus seems that 
there is no effect on lensing, contrary to the standard lore that light deflection 
is smaller in scalar-tensor theories. Actually there is a crucial difference since 
in scalar-tensor theory massive bodies do feel the scalar field. It follows that 
the gravitional constant measured in a Cavendish experiment today is not G* 
but Gn — G*j4o(l + a o) with a = din A/dip* and where a subscript indicates 
that the quantity is evaluated today. It follows that the dynamics of massive 
bodies, such as planetary orbits, determine GjyM and not G*MA 2 so that 

AG, MA 2 4G N M A9 GR 
M —be*— (l + a 2 )bc 2 ~ TTa 2 ~ ^ GR) 

as expected. Again, this shows that lensing alone cannot probe GR and that we 
need to compare different measurements. Note also that the gravitational con- 
stant (or more precisely the dimensionless number Gm 2 /he) varies with time 
so that extending this argument to the cosmological context is not straightfor- 
ward [57] . When the theoretical framework is specified then the post-Newtonian 
parameters can be computed (here 7 PPN = — 2a 2 / (I + a 2 ) and /3 PPN = 
a 2 /[2(l + a 2 ) 2 ]da/d^?* as long as the potential is such that the field is light 
on Solar system scales) so that the PPN constraints can be translated to con- 
straints on the parameters of the model. 

In conclusion, this simple extension of GR illustrates that we always have to 
introduce new fields in the theory so that we have to specify their nature (here 
a scalar field) and the ways they couple to the matter fields (here universally 
with the strength a). The distinction between a modification of GR and dark 
matter (or energy) is thus slight since in both cases we need to introduce new 
fields in our theory. The main difference lies in the fact that the amount of 
dark matter or dark energy is set by initial conditions (e.g. the amount of dark 
matter is fixed initially and determines the properties of the potential wells 
in which baryonic matter falls to form the structure or the amount of dark 
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energy is fixed by tuning some parameters and/or initial conditions so that it 
starts dominating today and the fact that p\ ■ Pcdm ■ Pb ~ 14 : 5 : 1 today 
calls for an explanation). In a modified GR model, the way standard matter 
generates potential wells or affect the dynamics of the universe is changed. 
Note however that the new degree of freedom are also gravitating so that in 
some models the distinction is even more subtle. 

Among the most studied alternative theories of gravity, let us mention 
scalar-tensor theories discussed above, f(R) theories which are, after a field 
redefinition, a sub-class of scalar-tensor theories, the DGP model [28] and the 
TeVeS [29] theory which is a relativistic version of the MOND [30] idea. 



2.4 Lessons for extending the tests to astrophysical scales 

GR is a well-defined theory of gravity with clear predictions so that the consis- 
tency of these predictions offers the possibility to test the theory in a model- 
independent way. This implies that we need various observables relating the 
same physical quantities (such as the mass in the example of the introduction). 

In cosmology, we can use almost the same observations as in the Solar 
system. Concerning light deflection, it cannot be measured (since the "unde- 
flected" position of the sources cannot be determined; but for the particular 
case of microlensing) and we will have to use the distortion of light bundles, 
that is strong and weak lensing. Also, and contrary to the Solar system, we can 
have access of the evolution of the energy of the photons, related to the time 
variations of the gravitational potential in the case of the integrated Sachs- 
Wolfe effect. The dynamics of massive bodies can be obtained from the large 
scale structure of the universe, which give an information of the growth of 
the structures and their velocity. Among the tests of the equivalence principle, 
only the test on the constancy of fundamental constants can be generalized. 

There are however limitations specific to cosmology. In particular, the cos- 
mological structures evolve with time and this contains an information on 
gravity but also on the properties of matter which are difficult to disantangle 
(for instance, our prediction on the shape of the galaxy power spectrum are 
different whether there exist massive neutrinos or not). This also means that 
we may have to take evolution effects into account. Also, cosmological data 
have to be interpreted in a statistical way so that we always have a depen- 
dence of the initial conditions that cannot be forgotten. Then, the description 
of the dynamics of the universe involves the Copernican principle so that the 
interpretation of our tests will depend on such a hypothesis. 

It follows that the tests that will be designed are indeed tests of GR but also 
depends on many other hypothesis so that they should probably be considered 
first as tests of the ylCDM model. 
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3 Galaxy and Cluster scales 

The first interesting systems for testing GR in astrophysics are galaxies and 
clusters. It is now well-established that, as long as one assumes GR to hold, 
their dynamics can only be understood by invoking the existence of dark mat- 
ter. 

The visible mass of spiral galaxies is rather concentrated so that Newtonian 
gravity predicts that the rotation curves should drop as r" 1 / 2 outside the 
bright part of these galaxies. But this has not been confirmed by more than a 
hundred rotation curve measurements [31j . Actually, in most spiral galaxies, 
and more particularly those with a high surface brightness, the rotation curves 
flatten at large distance from the center, — > const. Moreover, this velocity 
is correlated to the luminosity of the galaxy. This correlation, known as the 
Tully-Fisher law, states that the luminosity of the galaxy scales as , so that 
one expects that oc M, for the total stellar mass. This has provided the 
basis of the dark matter explanation: if the velocity is constant in the outer 
region of the galaxy, this means that the centripetal acceleration scales as 
a r oc r _1 and Newton's law implies that the gravitational potential scales as 
lnr. In the case of spherical symmetry, the Poisson equation implies that it 
should be sourced by a matter whose density profile scales as p(r) oc r~ 2 , as for 
an isothermal sphere model. Thus, each spiral galaxy must contain a spherical 
dark matter halo with a density profile scaling as r~ 2 at large distance. This 
reflects the discrepancy between two estimations of the mass: the luminous 
mass and the dynamical mass. 

To avoid such an hypothesis, Milgrom [3U] proposed a phenomenological 
modification, called MOND, that was able to account for the galaxy rotation 
curves [32], and more particularly to recover the Tully-Fisher law. MOND 
introduces a fundamental acceleration do, of the order of 1.2 x 10 _10 m • s~ 2 , 
such that the acceleration of any massive body is 



so that, at large distance, the Newtonian acceleration being GM/r 2 , the cen- 
tripetal acceleration is y/GMao/r. Since it is also given by v 2 /r, one deduces 
that v(r) -> (GMa ) 1/4 . In this regime, the gravitational potential behaves as 
\J GM ao In r instead of the standard Newtonian potential —GM/r. It follows 
that the deflection angle at large distance from the center of the galaxy is 



This value is the same as the one expected from GR, as long as one is in 
the halo. Indeed, if interpreted within GR, the presence of dark matter sug- 
gested by the rotation curves is confirmed by lensing observations. Therefore 
in MOND, an in any modification of GR, one must predict that a given mass 
generates a larger potential and a larger deflection angle than in GR. 

It follows that, one needs to estimate the mass of the galaxy, in a given 
theory of gravity, by different methods in order to check their compatibility. In 



a = ajVi if a > ao, 



a = ^/ajyao, if a < ap 



^MOND 



2wy/GMao 
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particular, different notions of mass needs to be distinguished [26]: the baryonic 
mass Mb assumed to be proportional to the luminous mass, or stellar mass 
M*; the total dynamical mass M t J™, estimated from the rotation curves; and 
the total tensing mass M^J t n determined by lensing observations. In the dark 
matter interpretation, and as well established by lensing observations, we have 

Dark Matter : M b < M t d y t n - M^ s . 

In particular, such mass estimates were performed in Ref. [33] using six strong 
lensing galaxies from the CASTLES database. The total mass was estimated 
from lensing while the stellar mass was estimated from a comparison of pho- 
tometry and stellar population synthesis. It demonstrates that dark matter is 
still needed (and that it is detected even in region where a > ao). In partic- 
ular this dark matter component cannot be explained by 2 eV neutrinos (see 
below) . 

If one assumes that the light deflection is given as in GR, then the previous 
equivalences told us that MOND predicts the same lensing as in GR within 
the dark matter halo. In particular the convergence at distance r is given by 
n(r) — r E /2r, where the Einstein radius is 

••^ 2 "(v) J ff < M0ND » "-*(t)'i£ (dm) - 

where the latter holds for a singular isothermal sphere with line-of-sight vel- 
city dispertion a v . While formally similar, these expressions have however an 
interesting difference [34] since a% scales as Mtot while scales as \fMl so 
that 

r E tx ^/Ml (MOND) r E oc M (DM). 

The scaling of the Einstein radius with the stellar mass was measured [34] 
using the RCS and SDSS surveys to show that r E ~ M°- 7i±0 - m . This seems 
in contradiction with the MOND prediction but the data used measurements 
of the shear at distances of some hundred of kpc, at which the environment 
effects can change the MOND prediction. It sets no constraint and the cold 
dark matter model since the fraction M*/M is not known. 

Indeed, MOND is a phenomenological description but not a field theory. 
As, we have seen earlier, the light deflection in scalar-tensor theories in smaller 
than in GR. This means that a MOND cannot derive from a simple scalar- 
tensor theory. It was realized (see Ref. 26 for more details) that this can 
be solved by coupling matter not to the metric g^ v but rather to a "physical 
metric" involving both a scalar and a vector field 

9 M v = e~ 2 V<> - 2U M U V sinli20. 

The first term is similar to what is performed in scalar-tensor theories and the 
second term, involving the vector field allows to reconcile light deflection 
with the GR prediction. This theory is known as the TeVeS (Tensor- Vector- 
Scalar) theory j 29 . Ref. [35] compares the TeVeS predictions to a large sample 
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of galaxy strong lenses from the CASTLES sample. Recently, Ref. [36] com- 
pares the predictions of TeVeS for both galaxy rotation curves and strong lens- 
ing (for high and low surface brightness galaxies) concluding that TeVes, in its 
simplest form, cannot reproduce these data consistenly. The analysis [37] of 
galaxy-scale strong lensing from the Sloan ACS (SLACS) survey indicates that 
I^ppn _ i| < 5 x 10 -2 . However this work emphasizes that setting constraints 
with such systems requires to know the properties of the lensing galaxies with 
a great accuracy, much greater than at present. The comparison of the stellar 
velocity dispertion to measurements of the Einstein radius allows to constrain 
7 PPN , reaching [35] 7 PPN = 0.88 ±0.05 on kiloparsec scales at 68% C.L. while 
an early analysis based on 14 systems only gave [39] 7 PPN = 0.93 ± 0.1. 

On cluster scales, various estimates of the mass can be obtained by lensing 
(strong and weak), X-ray emission that characterizes the intracluster (bary- 
onic) medium, and the SZ effect which gives an information of the electron 
distribution. By comparing these distributions, one can compare the location 
of the gas and the gravitational iso-potential. Earlier analysis used the com- 
parison of X-ray and strong lensing [40] and then weak lensing [41j leading to 
the conclusion that the Poisson equation should be valid, within a factor 2, up 
to scales of 2 Mpc [42] . 

Interesting conclusions arise from the study of the colliding galaxy clusters 
1E0657-56 (z = 0.296). In this system a smaller cluster, known as the "bullet 
cluster" , has crashed through a larger one and their intracluster gas has been 
stripped by the collision. On one hand, weak lensing shows that the lensing 
mass is concentrated in the two regions containing the galaxies rather than 
in the stripes containing the baryonic matter 43J. A similar observation [44] 
was made with the merging galaxy cluster MACS J0025. 4-1222 (z — 0.586) for 
which the emitting gas, traced from its X-emission, is clearly displaced from the 
distribution of galaxies (from lensing). The rich cluster Abell 520 (z = 0.201) 
also exhibits the same properties [IS] and contains a massive dark core, as 
deduced from lensing mass reconstruction, that coincides with the central X- 
ray emission peak. The analysis [46] of the cluster Abel 478 demonstrates that 
the X-ray, SZ and weak lensing data perfectly agree with a dark matter model 
(but does not prove they cannot be reproduced by a MOND model). 

This seems to be a proof of the need of dark matter since, being collision- 
less, it continues to be located around the bullet, contrary to the baryonic 
gas. This was confirmed [47j by the reconstruction of the mass distribution 
from both strong and weak lensing. However, it seems that MOND could 
accomodate these observations in particular because the original TeVeS ver- 
sions make different prediction when the system is not spherically symmetric. 
Ref. [48] showed that it was possible to design a multi-centred baryonic system 
eproducing the weak-lensing signal of 1E0657-56 with a bullcr-likc light dis- 
tribution. The same authors [35] then realized that a purely baryonic MOND 
model cannot accomodate the data and that the bullet cluster was dominated 
by dark matter whether one uses GR or MOND. In the latter case, it would 
require massive neutrinos with m v = 2 eV. 
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The existence of such neutrinos seems however problematic. From the study 
weak lensing for 3 Abell clusters and 42 SDSS clusters, it was concluded that 
MOND cannot explain the data unless some dark matter is added and this 
dark matter cannot be accounted for by massive neutrinos [50] . This was con- 
firmed [51] by the confrontation of strong and weak lensing from the HST 
Wide-Field Camera, excluding the dark matter to be neutrinos with mass in 
the range 2-7 eV. 

In conclusion, it seems that MOND and TeVeS have difficulties to repro- 
duce the observations of the distribution of dynamical, baryonic and lensing 
masses. Indeed none of the above mentioned results demonstrate that MOND 
is ruled out. They are analysis that show that the data can be consistently 
interpreted assuming GR and the existence and dark matter. One of the main 
difficulty to use these observations as a direct test of GR is the complex ge- 
ometry of the systems that are used. 

4 Cosmological scales 

The construction of a cosmological model relies on the choice of the theory 
of gravity as well as on our understanding of the fundamental interactions 
of nature. However it also involves other hypothesis, such as the Copernican 
principle which states that we are not seating in a priviledged place of space. 
Under such an hypothesis, and whatever the theory of gravity, the universe on 
large scales can be described by a Fricdmann-Lemaitre spacetime with metric 

ds 2 = -dt 2 + a 2 (t) 7 yda;Mx J , (5) 

where t is the cosmic time, a the scale factor and 7^ the spatial metric on 
constant t hypersurfaces. If one assumes that GR is a good description of 
gravity then the dynamical equations of such a model derives from the Einstein 
and the conservation equations 

= SnGT^, V^T"" = 0, (6) 

where G M „ is the Einstein tenor and the total stress-energy tensor. Among 
this class of models our reference model is the vlCDM model which includes a 
cosmological constant and cold dark matter and assume that initial conditions 
are consistent with the prediction of slow-roll inflation. Such a model is in very 
good agreement with most of the astrophysical data and it is self-consistent. 
But the fact that the dark sector represents 95% of the matter content of the 
universe and the cosmological constant problem drive us to test the hypothesis 
of our model, and in the first place the Copernican principle and GR. It is 
indeed difficult to anbandon these two hypothesis at the same time so that all 
the studies aiming at testing GR in that context assume that the spacetime 
metric remains of the form Also most of them still include dark matter 
and aim at replacing dark energy by a modification of GR. 

Two roads can be followed. Either one defines a class of gravity models that 
contains GR in some limit and then confronts it to cosmological data to see 



12 



how close from GR, in this particular space of theories, the theory of gravity 
should seat. As an example, this was performed in depth for scalar-tensor 
theories (|4j) for which the implications of the background dynamics [52], of 
the cosmic microwave background [35], primordial nucleosynthesis [S3], weak 
lensing [27], and local constraints [54], even allowing for extensions to a non- 
universal coupling of dark matter [55] , were all studied. Or one tries to quantify 
the allowed deviations from the reference model while being as much as can 
be model- independent. The strategy is then to exhibit consistency relations, 
analogous to Eq. ([T]), betwen different observables, which must hold in our 
ylCDM reference model, as first proposed on the particular case of the Poisson 
equation in Ref . 3 J . As explained in Ref . [56] , the modification of our reference 
framework can be classified in universaltiy classes who have specific signatures 
and different tests can favour or disfavour some classes of modification. 



H 2 = — (p + Pde) - — , - = — (p + 3P + p de + 3P de ), (8) 



4.1 Background dynamics 

The dynamics of the background spacetime is dictated by the Friedmann equa- 
tions 

where H = a/a is the Hubble function and K = 0, ±1 is the curvature of the 
spatial sections. From the last equation, the recent acceleration of the cosmic 
expansion implies that (p + 2>P) < if GR is a good description of gravity. 

At the background level, a modification of GR or the introduction of a 
dark energy component instead of the cosmological constant will change the 
Friedmann equation and can be taken into account in a effective way simply 
in terms on an effective fluid 

8ttG. , K a AttG 
— (p + p dc )-^, - = 

O Li Li O 

where p dc and P dc can depend on H and its derivatives, as e.g. for scalar-tensor 
theories or DGP. This allows to define the equation of state of the dark energy 
from H (see Ref. [54]). 

Indeed, without an explicit model, the extra-terms p de and P de are not 
known. Besides we known that, since they arise from the existence of a new 
degree of freedom, there must exist an associated equation of evolution [57] , 
The standard approach is to postulate that they are related by an equation of 
state, pdc = w dc P dci and the most commonly used ansatz is [58,59 

w de (z) = w + — — w a , (9) 
1 + z 

with wq and w a constant. It is thus clear that the background dynamics cannot 
distinguish between a modification of GR and a properly tuned dark energy 
model. This lies in the fact that the only quantity at hand is H (z) and most 
of the models of the literature can be tuned to reproduce the same function 
(see Ref. [S3] for explicit examples). No null test of GR can be constructed 
with background data since they all are functions of H(z). 
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4.2 Linear perturbation theory 
4.2.1 Standard ACDM 

As long as one assumes GR to hold and consider an almost Fricdmann- 
Lemaitre spacetime, the evolution of perturbations is well understood, see e.g. 
Ref. [60]. On sub-Hubble scales, focusing only on scalar perturbations which 
are dominant at late time, the space-time metric can be written as 

ds 2 = -(1 + 2<P)At 2 + (1 - 2y)a 2 {t) ltJ &x % &x\ (10) 

where <P and <P are the two gravitational potentials. 

The evolution equations on Hubble scales are given by the conservation of 
the matter stress-energy tensor (continuity and Euler equations) 

6 m = - — , ni + H6 ni = --A<P, (11) 
a a 

which leads to the standard second order evolution equation 

5 m + 2HS m - = 0. (12) 
Among the Einstein equations, we can keep only the Poisson equation 

A& = AirGp nl a 2 S m (13) 

and 

<5-<f = (14) 

that arises from the fact that the matter anisotropic stress is negligible. This 
shows that the two gravitational potentials have to coincide, which is related 
to the fact that 7 PPN = in GR, their spectrum has to be proportional to 
the matter power spectrum, as firt pointed out in Ref. 3 and Eq. (|11[) implies 
that 

e m = -fS m , with / = ^r^, (15) 

dm a 

where we have decomposed the density contrast as <5 m = D(t)e{x) where e 
encodes the initial conditions. It was shown [61] that, for a flat CDM model, 

/ ~ (16) 

was a good fit. The important feature here is that if GR holds, then the growth 
of structures is completely determined by H(z). We forsee that one can check 
the compatibility of background data (such a distance-redshift relations, e.g. 
from SNIa) and large scale structure data, whatever the parameters entering 
the equation of state © or any other parameterisation (see below). 
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Jf..2.2 Extension 

In order to construct tests of GR in this regime, one needs to construct the 
most general extension of this set of evolution equations, still assuming we are 
dealing with a metric theory of gravity. 

In the case of dark energy alone, one needs to consider the effect of its 
stress-energy tensor, which can have non-vanishing anisotropic stress and den- 
sity contrast, contrary to a pure cosmological constant but the equations of 
evolution of the other fluids are not modified, since otherwise this new com- 
ponent would be coupled to the standard matter non minimally. In GR mod- 
ifications, there exists a new long range force and the evolution equation of 
matter will be of the form V^T/ 1 " = f", where is a force term between 
the standard matter fields and the new degree of freedom. The way such force 
term appears in the equation and its relation to the Einstein equation is not 
obvious to describe in full generality while being model- independent. As an 
example, consider scalar-tensor theories of gravity. In the Jordan frame, the 
equations of motion of the standard matter fields are not modified so that 
ff 1 = but, performing a conformal transformation, the same theory, written 
in the Einstein frame, involves a force f % = a^^Tid^ip* that will appear even 
at the background level in the continuity equation. In this particular case, it is 
well understood that the modification of gravity appears as a time-dependent 
modification of the Newton constant in the Jordan frame while it is seen as a 
universal time-dependent modification of masses in the Einstein frame. Indeed, 
if the new force is not universal, it probably involves that mass ratios will be 
time-dependent, which can be tested [6]. Thus, we assume that we are working 
in the equivalent of the Jordan frame so that we assume that /P = and that 
there is no creation of matter. The spatial component of the force can how- 
ever be non vanishing and enters the Euler equation. Let us stress that while 
important from a physical point of view [62,63 , this is not dramatic from a 
phenomenological point of view since only the source term of Eq. (|18|) below 
will be changed. On the other hand, we shall have an equation of evolution 
for the new degree of freedom that shall also have a source term proportional 
to the matter stress-energy tensor [62]. Unfortunately, this equation remains 
unkown until we specify the model. 

As long as we stick to linear perturbations, these extensions can be imple- 
mented by modifying the previous equations as 




m + H9 m = --A<P + S dc 



a 



(17) 



which leads to the standard second order evolution equation 



'L + 2HS m - -^A<P = S dc . 



(18) 



a- 



The term S de (k, a) encodes the new long-range force between the new degree 
of freedom and the standar matter (and dark matter!). 
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Then, we need to write down the Einstein equations. First, we can gener- 
alize the Poisson equation, written in Fourier space, as 

- fc 2 <f = 47rGF(fc, H)p m a 2 S m + A de . (19) 

The first term F(k, H) accounts for a scale dependence of the gravitational 
interaction while Zide accounts for a possible clustering of the new degree of 
freedom, and in particular of dark energy if there is no modification of GR 
(this shows at this stage, that the Poisson equation can be modified without 
modification of GR if dark energy can cluster; also care needs to be taken in 
the case of massive neutrinos which can enter on the r.h.s. of this equation, 
see e.g. [53]; Ref. proposed a interesting example of a clustering dark 

energy model mimicking the DGP model). Then, there is the possibility to 
have an effective anisotropic stress so that 

A($ -9) = Tide- (20) 

It follows that the deviation from GR is encoded in the four functions (Sd e , F, 
Ade, Ude) which, in the case of the /1CDM model, reduces to (0, 1, 0, 0) and, in 
the case of dark energy to (0, 1, A do , Tide), even though in most cases Zide and 
Tide are negligible. Their expression for quintessence, scalar-tensor and DGP 
models can be found in Ref. [56] . For the same reason that, in the case of a 
GR modification, at the background level Pde and pde can depend on H, Sde, 
Ade and Tide can depend on <P and while F is a function of the background 
quantities only. 

Equations (|17I19[) imply that the matter density evolves as 

L + 2H5 m - 4nGF(k, a)p m a 2 6 m = C dc (21) 

with Cdc = (^dc + I7dc)/a 2 — Sdc/a so that the matter power spectrum is 
expected to be deformed in shape, mainly because of the fc-dependence arising 
from F and from the source term Cdc- 

This idea to construct such a post-vlCDM parameterisation on sub-Hubble 
scales was first proposed in Ref. [56] . following the analysis of the particu- 
lar case of scalar-tensor theory by Ref. [57] • Several similar approaches were 
then designed in Refs [B^rBTVB^IJBWT^rTTirTWTWTW^r^ which went further 
in determining the relations with observational data (see discussion below) . In 
particular Ref. [67] uses, instead of Tide, a parameter rj defined by 

<P = (1 + ri)$f, 

so that it is a generalisation of the post-Newtonian parameter 7 PPN , as first 
proposed in the particular case of scalar-tensor theories in Ref. [27] . {$> ^ # 
in the Solar system is an indication of a modification of GR because one is 
dealing with vacuum solution of Einstein equation; again, this is not the case in 
cosmology since dark energy can have an anisotropic stress). Instead of using 
the the set (/, A de ), most analysis, including Refs. 66 , 671 17511731175] . assume 
that the Poisson equation is modified to 



-k 2 <F = 4irGQ(k, a)p m a 2 S, 
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involving only one new function Q. In most cases at hand, this is a good ap- 
proximation, but in full generality we shoud distinguish the large-scale modi- 
fication of GR and the clustering of possible new degrees of freedom (see e.g. 
Ref. 65 , 66J ) . In particular A^ e could be an independent random variable, not 
proportional to <5 m . The limit zido = corresponds to pure modification of 
GR with negligible effect of the new degree to the total stress-energy tensor 
while F = l,Z\dc 7^ corresponds to models of clustering dark energy (and 
may also incorporate the effect of massive neutrinos; see e.g. [77] for a study 
of this degeneracy). 

As such, this description is not complete since we have no equation to de- 
scribe the evolution of the new degrees of freedom. At the background level, 
this gap is often filled by assuming a parameterisation of the dark energy 
equation of state. As we shall see, two roads may be followed from this point: 
either one parameterized the unkwon functions that appear here or one con- 
struct null-tests. 



4.3 Cosmological data 

The previous analysis shows that in order to constrain the deviation from GR 
with the large scale structure of the universe we need to be able to extract 
information on the distribution of the four variables (<P, SP, <5 m , 9 m ) which are 
not directly observable. Let us summarize briefly some of the observations that 
turn to be useful. We call Pxy(k, z) the 3-dimensional power spectrum of the 
fields X and X at redshift z (or equivalently time t or distance x) defined 
by (A(k,z)F(k',z)*) = (27r) 3 F X y(fc,z)(5^(k - k') and C XY (e,z) their 2- 
dimensional (or angular) power spectra. 

4-3.1 Background data 

Background data usually include the luminosity distance-redshift relation pro- 
bed by SNIa which provides a handle on H(z) up to z ~ 1.5, the angular 
diameter distance mainly from the tangential component of the B AO measured 
by their imprint on galaxy distribution. 

4-3.2 Large scale structure 

The clustering of galaxies is one of the oldest measures of the properties of the 
large scale structure. The galaxy power spectrum P gg is the simplest measure of 
the correlations in the galaxy number density n g . In general the distribution of 
galaxies is biased with respect to the mass distribution and it is often assumed 
that they can be related by 

8 g ^ = M m + (22) 

Tin £ 
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where b\ and 62 are biais parameters. It is expected that the biais is a com- 
plicated function of time and of the masses of the halo hosting the galax- 
ies [78], since it encodes in some way all the process of galaxy formation. 
We shall assume here, for simplicity, that we restrict to a linear biais so that 
P gg {k, z) = bi(z)Ps JXL s m (k, z). Imaging data with photometric redshift provides 
a measurement of the angular power spectrum of galaxies, which is a simple 
projection of the three-dimensional power spectrum 



c ag (t) 



w 2 ( x ) 



Sk(x) 



p 

2 99 



Sk(xY 



(23) 



Sk being the comoving angular distance and where W g is the normalized 
redshift distribution of galaxies in the sample. 

However, the redshift-spacc position of any galaxy differs from its real space 
position due to its peculiar velocity. The density contrast in redshift space, 5 g , 
and in real space, 5 g , can be related by imposing mass conservation |79j . In 



the linear regime, this leads to Si 



/i 8g/aH, where /i is the cosine of 



the line-of-sight angle so that the redshift-space power spectrum is 



P s gg {k, i x) = P gg {k)+2^ j P g e g {k) 



i 2 H 2 



Pf> B eM. 



(24) 



It is thus commonly modelled as [751179 



P 9a {k) + 2^P geg {k) 



^Pee 



(k) 



F 



H 2 {z) 



(25) 



where F is a smoothing function and a v is the 1-dimensional velocity disper- 
sion. The angular dependence enables to separate the different components [5U] 
to get a measurement of the three spectra P gg {k, z), P g g g (k, z) and Pg g e g {k, z), 
in particular from the SDSS [5T] and 2dF [55] galaxy redshift surveys. These 
low redshift analysis were extended to z ~ 1 in Ref. [53] using the VIMOS-VLT 
Deep Survey (VVDS) [84] . 

Indeed, in the standard ylCDM, and in the linear regime, we have that 
g = —aS m so that its growth rate is Dg g = —aD so that 9 g = aHf8 m . In that 
limit P g = aH/3Pgg and Pe g e g — a 2 H 2 p 2 P gg with j3 = f/b, and the three 



spectra are not independent. 



4-3.3 Weak lensing 

Gravitational lensing offers various posibilities. As previously, we restrict our 
analysis to metric theories of gravity. 

First, either in the strong or weak regime, it can probe the sum of the two 
Bardeen potentials, # + <P\ Weak lensing surveys use the observed ellipticities 
of backgroung galaxies (and more particularly the correlation of their shapes) 
to reconstruct a map of the cosmic shear, which can then be used to determine 
the convergence n [85 . As long as photons travel on null geodesies and the 
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w(x) = s K (x) I w g ( x ') J * U ,„ A X) d x >. 



geodesic deviation equation holds, the distortion of the shape of background 
galaxies can be computed from the Sachs equation [86] leading to 

Xi) = \j W(x, xi)M* + *)Ax (26) 

for sources located in a bin centered round a redshift z% and with x% = xi z i) 
with 

w{ X , Xi ) = SK{x) " K ( {x :- x) . 

OK{Xi) 

The convergence power spectrum for two sets of galaxies centered around Zi 
and Zj, as can be obtained by a tomographic survey, is thus 

C KK (£,z l:Zj ) = W( X ,Xi)W(x,Xi)P$+# (k = s^ftj'X) dx - ( 2? ) 

Until we have data allowing for the use of tomography, we have only access to 
the shear power spectrum averaged on the source redshift distribution, P KK (£) 
which is given by the same expression but with the window function 

^ S K (x' - X) 
Sk(x') 

In conclusion, this allows to constrain the power spectrum P$ + &(k, z). Note 
that the analysis of weak lensing requires in fact to know the non-linear power 
spectrum but the latest data from the CFHTLS [87] reach large angular scales 
(9 > 30 arcmin.), which allows to work in the (quasi) linear regime, where 
theoretical predictions for the time evolution of the power spectrum are more 
reliable. Note also that the convergence power spectrum is often expressed in 
terms of the matter power spectrum, making use both of the Friedmann and 
Poisson equations. Indeed, the goal here is to relate the observables to their 
primary perturbation variables without using any equations. 

Second, one can use galaxy-galaxy lensing, which arises when the deflect- 
ing and source galaxies are aligned, giving rise to a mean tangential shear 
around foreground galaxies. This will thus give an information on the corre- 
lation between the galaxy distribution and <P + \P through the angular power 
spectrum 

C gK (i,Zi,Zj) = l - J W g ( X ,Xi)W( X ,Xj)P g ,*+* (k = s^y.x) dx- (28) 

This was for instance measured from the SDSS galaxy survey. Note that the 
magnification biais [78] can also help to extract some correlations since 

S g = bS m + (5s - 2)k 

where s — dlogTV/ dm is the logarythmic slope of the number count-magnitude 
function. This induces distortions 88,89] that can also be used to test GR. 
Note that it was also proposed that the weak lensing of standard candles 
(SNIa, or GW sirens) can be used to measure the cross correlation between the 
magnification \i and S g [9"0"ll9"T] . In particular C MM and C 9 ^ contain informations 
similar to C KK and C gti respectively. This has not been investigated yet. 
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4-3.4 Integrated Sachs- Wolfe effect 

The observation of the cosmic microwave background temperature anisotropics 
gives numerous important informations for our cosmological model, among 
which the initial power spectrum for the perturbations. While propagating 
from the last scattering surface to us, the energy of the photons changes due 
to the fact that they cross structures in formation, and thus propagate in a 
spacetime where <P and & are not constant. This induces a direction-dependent 
temperature change, known as the integrated Sachs- Wolfe effect [50] 

AT f ■ 

— = J {$ + &)a( X )dx, 

where the integral is performed along the photon geodesic. This ISW effect is 
correlated with the galaxy distribution with angular power spectrum 

C g ,isw(t) = J P g>[ S^) (k = J^yx) ^d x . (29) 

This has been detected [92] by cross-correlating the CMB anisotropies to 
galaxy maps. 



4-3.5 Conclusions 

We see that astrophysical observations allow to measure many correlations 
between the perturbation variables. In order to construct tests of GR, one 
needs to relate 5 g to <5 m and thus understand the biais (and most importantly 
contrain its scale dependence) . The example given above are the most promis- 
ing to implement the tests of GR but, indeed, there exist many other ways to 
measure these quantities. 

We also need to keep in mind that each of this method has its own sys- 
tematics and limits. We cannot discuss this issue here, but it is central when 
actually deriving constraints. 



4.4 Growth of matter perturbations 

The first effect of a modification of GR is to change the growth of density 
perturbation. In the /1CDM model, Eqs. (|lltil2[) imply that the growth rate 
D evolves as 

D + 2HD - iirGp m D = 0. (30) 
This equation can be recast in terms of a as time variable [60 as 

D» + (^ + 3 -)d>= 3 -^D, (31) 
V da a J 2 a 
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from which it can be checked that D = H is a solution so that the growing 
mode can be obtained as 

D = l^-n m0 r dM (32) 

which implies that if H(z) is known from background observations, such as 
SNIa, then D(z) is fixed. There is a rigidity between the expansion history 
and the growth rate. 

The growth rate can be parameterized phenomenologically as [931194] 

Then, if GR is not modified, the index 7 can be computed once H(z), or 
equivalently the dark energy equation of state, is known and it was shown |95[ 
[96] that 

7 = 0.55 + 0.05[1 + w dc (z = 1)]. 

While being a good test of dark-energy model with a smooth energy distribu- 
tion, it is not clear whether it can be considered as a test of GR. In particular, 
we can imagine that dark energy has an anisotropic stress or is clustering while 
GR is not modified (i.e. F = 1, Sd c = 0, 77d c 7^ 0, Z\d c 7^ so that <Sdc 7^ 0) 
so that one accomodate a value of 7 by some properly designed model. To 
finish, such a parameterisation is too restrictive since it does not include the 
scale-dependence that is expected from the modification of GR (see however 
Ref. [76]). The extensions to include the super-Hubble regime were considered 
in Refs. 69,97,76 . It was recently proposed [HE] that galaxy cluster velocities, 
measured from the kinetic SZ effect, may allow for a measurement of 7. 

Several studies concentrate on a pure modification of the Poisson equation 
so that only the term F(k,H) is modified. The effect of such a fc-dependent 
term on the power spectrum was first studied in Ref. [3] where the function F 
was assumed to reproduce the effect of higher-dimensional gravity, as described 
at the time. In that particular case, where the only parameter is the length 
scale r s , it was shown [100] that the cosmic shear [99] 3-point function implies 
that r s > 2/i _1 Mpc. 

Similar analysis in the case of a Yukawa type modification of the gravita- 
tional potential, 



<£(r) =-G J \ 



P(r') 



1 + a 1 - e 



,-|r-r'|/A 



were then performed. With such a potential, the function F(k, a) entering the 
Poisson equation is given by [1011ll02j 

Such a modification causes the rate of growth to depend on k so that the scale A 
shall have an imprint on the power spectrum (see also Ref. [103j for a general 
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argument on the shape dependence). Interestingly, assuming an Einstein-de 
Sitter bacground cosmology, Eq. (|2l|) can be solved analytically in terms of 
hypergeometric function [101(1102] to give the growing mode 

5 - V25 + 24a 5 + V25 + 24a 9 2 \ 
8 ' 8 '4 5_S J 

with s = a/kX. 

Ref. [101] considered the case of a Einstein-de Sitter background and ana- 
lyzed the SDSS and 2dFGRS data up to k ~ 0.15/i/Mpc, leading respectively 
to the constraints a = 0.025 ± 1.7 and a = —0.35 ± 0.9 at a la level. A sim- 
ilar analysis was performed in Ref. [102] who used the Peacock and Dodds 
procedure [104j to describe the non-linear power spectrum. The analysis of 
the SDSS data sets the constraints —0.5 < a < 0.6 (resp. —0.8 < a < 0.9) 
for A = bh- 1 Mpc (resp. A = lOftr 1 Mpc). The analysis was extended in 
Ref. [105] by performing both second order perturbations and N-body simula- 
tions to construct a mock galaxy catalog. Ref. |106] extended these analysis by 
allowing a modified expansion rate, which should be the case if GR is modified. 
They also showed that the modification of the shape of the power spectrum 
is almost degenerate with the effect of massive neutrinos. Notice that a com- 
bined analysis [107] using CFHTLS weak lensing data and the SDSS matter 
power spectrum estimated from lumianous red galaxies found no sign of devi- 
ation from GR on scales ranging between 0.04 and 10 Mpc. Even though this 
analysis used both matter distribution and weak lensing, it only constrained 
the shape of the power spectrum without implementing the consistency chek 
proposed in Ref. [3]- 

N-body simulations with such a Yukawa modification of the gravitational 
potential were performed in Ref. [108] who concluded that the gravitational 
evolution is almost universal, at least for A in the 1 — 20 Mpc range so that 
the Peacock and Dodds approach [104j can be adapted to get an analytical fit. 
It was extended by the simulations of Ref. [109] which include the possibility 
of an anisotropic stress and considered the case of DGP models with r s = 
(5, 10, 20)h~ 1 Mpc. To finish, the spherical collapse model and the estimate 
of the abundance of virialized objects was considered in Refs. [11011111 . The 
scale dependence of the growth rate was proposed [112j to be studied in terms 
of 

which remains close to for any smooth dark energy model. In particular it 
can be measured from future redshift surveys. 

These studies allows to understand the effect of the modification of the 
Poisson equation, which is expected to be generic in any deviations from GR, 
and absent in all models of pure dark energy. They are thus very instructive but 
note that the background cosmology is in general not modified in a consistent 
way. 



5 + (k,a) = s 2 Fi 
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4.5 Testing GR on cosmological scales 

There have been two main approaches to using cosmological data to constrain 
deviations from GR. 

4-5.1 Parameterizing our ignorance 

In the first approach, one tries to use the generalized set of perturbation equa- 
tions (|17H20p in order to compute various cosmological observables and com- 
pare them to astrophysical data. The main problem, as mentioned earlier, is 
that such a parameterisation cannot be complete unless the physics of the 
new degrees of freedom is know. We thus have two possibilities. Either one 
compute explicitely these terms in some classes of theories such as f(R), DGP 
or scalar-tensor [56 l l67 l l68] or one specifies some ansatze for these functions, in 
the same spirit as we introduced the parameterisation ([9|) for the dark energy 
equation of state. 

For instance, Refs. [671173] assume that the function S = Q(a,k)(l + 
T)(a,k)/2) can be expanded as U = 1 + E^a with Sq constant. The effect 
of the modification of GR is taken into account through a parameterisation 
of the form (|33|) , with 7 constant so that one ends up with 4 constant extra- 
parameters (wo, w a , 7, £0) besides the standard cosmological parameters, the 
ACDM model corresponding to (wq, w a , 7, -Eb) = (—1,0,0.55,0). Such a pa- 
rameterisation was then used to discuss the sensitivity of various probes. 
Clearly, this choice misses a possible scale-dependence of F (or Q) which is 
generically expected if GR is modified |103j . This issue was recently adressed 
in Ref. |113j which proposes to expand the two unknwon functions Q and rj as 

X(a, k) ~ X (a) + Xx{a)aH/k. 

Ref. [114] chooses the same two functions as functions of k and a and proposes 
different ansatze for their functional form in order to study the potential of up- 
coming and future tomographic surveys to constrain them. Ref. [68] proposed 
a parameterisation that depends on the scale. On the other side Refs. [7T|l72] 
focus only on the function r)(z), that they call w{z), in order to infer its influ- 
ence on CMB anisotropy spectrum and weak lensing. zu was chosen to scale as 
^oPde{z) / p m (z), assuming a /1CDM evolution for the background cosmology. 

Note that in the standard ACDM, we must have that ^(k, a) = ^(k, a) — 
—3Q m (a)(Ha/k) 2 S m (k,a)/2 and m (k, a) = — fS m (k, a). Thus, instead of pa- 
rameterizing the unknown terms that enter the perturbation equations, we 
may think to parameterize directly their solution, for instance, as 

#(k, a) = -^/2 m (a)(iTa/fc) 2 £ m (k, a) [1 + c*(ft, a)] , 
!F(k, a) = -^f2 m (a){Ha/k) 2 6 m {k, a) [1 + c»(fc, a)] , 

and 

m (k, a) = -f5 m (k, a) [1 + ce(k, a)) , 
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where <5 m (k, a) is supposed to scale as <5 m (k, a) = S m (k, a,i n )D(a) [1 + c$(k, a)] 
and then finding physically motivated ansatze for the functions Cj. Such an 
alternative parameterisation was considered in Ref. [70j . 

In these approaches, the game is thus to replace free unkown functions by 
a set of parameters in order to be able to compute the different cosmologi- 
cal observable and then compare them to data. This allows in particular to 
understand the accuracy with which they can be constrained by forthcoming 
experiments, the main difficulty being to find the most relevant set of parame- 
ters that reproduce a large class of theories. Ref. [73] utilised the large angular 
scales (9 > 30 arcmin.) weak lensing data from the CFHTLS (in order to work 
in the linear regime), BAO and SNIa data to get constraints on (J7 m o, Uo,l) 
consistent (e.g. related to the parametrisation of the equation of state) with 
the standard ylCDM. Note also that it is important that the background and 
perturbation dynamics be consistent since they derive from the same modifi- 
cation of GR, an issue often overlooked. 

4-5.2 The art of correlating 

A probably better idea to obtain constraints on deviation from GR is to con- 
struct null tests. Such tests are based on the simple fact that once the theory is 
completely specified, there must exist consistency relations, between different 
observables, in a similar way as the Solar system example of the introduction 
led to the consistency relation (TTJ). Indeed, any departure from such a relation 
would indicate that some hypothesis of our model are not correct and that the 
theory needs to be extended, without telling how. Such tests are null tests, in 
the spirit of "traditional" physics in which a reference model is confronted to 
observations in order to determine the limits of its validity. 

The use of cosmological data to perform such tests was first proposed in 
Ref. [3] who focused on the Poisson equation (fT3|) . If such an equation holds 
then the power spectra of the gravitational potential Pj(fc) and of the matter 
distribution P$ m (k) must be related by 

fe 4 Pp(M) = ^f2 m0 H^a- 2 P 5m (k,a), 

whatever the cosmological scenario. This means that the scale dependence of 
the two pectra are related in a very specific way. In particular, if the Poisson 
equation is modified, the change of the shape of the matter power spectrum 
is, as we saw on the example of a Yukawa potential above, model dependent 
but the fact that the two spectra differ is a model independent conclusion. In 
particular such a relation can be tested by comparing weak lensing data to 
galaxy survey, if the scale dependence of the biais is mild, as expected from 
numerical simulation, since C KK and P gg give access to P& and P,5 m . Note also 
that, it has a trivial generalisation if the fields are all proportional to the same 
stochastic variable (which is the case for adiabatic initial conditions) then 
P<M m = \/P^PsZ, which again can be tested using galaxy-galaxy lensing. 
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Similar rigidities were exhibited between the background dynamics and the 
growth of the large scale structure. For instance, the growth equation (|3ip. 
valid for a vlCDM, and thus when GR hold, can be recast |115pil6j (see also 
Refs. |1171lll8j ) as a first order equation for H so that H(z) can be inferred 
from background data and perturbation data independently, or cquivalcntly 
the equation of state of dark energy ([9]) and the parameter 7 defined in 
Eq. (f33|) are not independent when the dark energy is assumed to remain 
smoothly distributed. This was implemented in the analysis of Refs. |119pi20] 
who introduce two dark energy equations of state, one for the evolution of the 
background geometry and the other governing the growth. Using SNLS-SNIa, 
2dF and SDSS galaxy redshift survey, CMB data and CTIO-lensing survey, 
they concluded that the two determinations of Qa were consistent and that 
the two constant dark energy equations of state have also to agree. These anal- 
ysis consider only the effect of the growth factor and no other modification 
is considered. Another implementation performs a model-independent recon- 
struction of the growth rate from distance measurements and then compares 
to growth measurements |12ipi22j . Ref. |123j proposed a similar consistency 
test of the /1CDM using low and high redshift SNIa survey by estimating fl m 
in three different ways (background geometry, growth, and shape of the power 
spectrum), all agreeing with the canonical value 0.25. 

The original idea of Ref. [3] was extended to multiple cosmological probes. 
Ref. [124) proposed to use the galaxy-velocity correlation and the galaxy- 
galaxy lensing, which give access to (S g 9 m ) cx bf(S^ n ) and (6 g n) cx b(5 m A(<l> + 
!?)) so that the ratio of these two quantities is expected to be independent 
of the biais, at least in the regime of linear biasing. An estimator, Eq based 
on the ratio of these two quantities, was constructed and it was demonstrated 
that it can distinguish a large class of models. Ref. [66] worked out the re- 
lations between the various observables, including a discussion of quasilinear 
effects. It was also shown that the clustering of dark energy can mimic features 
of a modification of GR and investigated the way to combine data in order 
to distinguish the two effetcs. Refs. [731175] designed consistency checks based 
on the redshift-space power spectrum and weak lensing in order to constraint 
the ratio and the Poisson equation. Ref. [125] proposed an estimator to 
measure the ratio of the two gravitional potentials, again using weak lensing 
and redshift-space power spectrum. Ref. |126j proposed a method to extract 
the effect of a modified Poisson equation and Ref. [113] analyzed the combi- 
nation of imaging and spectroscopic surveys. Ref. [127] proposed to use the 
ISW-structure correlation to constrain the growth rate of the density, and in 
particular its scale-dependence. 

All these works are thus starting from the constitutive relations that exist in 
a ylCDM, and thus assuming GR valid to construct from large scale structure 
survey some tests that will indicate the violation of one of these relations. 
They often construct estimators that are probed by using some extensions of 
GR (DGP, f(R), scalar-tensor) and in oder to forecast the power of coming 
surveys to distinguish between them. 



25 



5 Conclusions 

This review has presented the tests of GR on astrophysical scales, but more 
generally of the /1CDM model, based on the large scale structure and the 
global dynamics of the universe. In particular, it is important to make the 
distinction since some supposedly tests of GR proposed in the literature are 
in fact only tests of the /1CDM model. Maybe the first answer these tests will 
give is whether there is a need for new physical degrees of freedom in our 
model and then start to caracterize the nature and the couplings of this field 
with standard matter (and also dark matter). 

As explained, future surveys will allow to map weak lensing, galaxy distri- 
bution and velocity on sub-Hubble scales with high accuracy. This will allow 
to construct many consistency checks of our cosmological model, and in par- 
ticular of GR. A multi-probe approach will allow to have a better control of 
systematics which affect each probe. 

Today, data shows no deviation from the ylCDM model, and thus from 
GR on large scales. On galactic scales, the debate between dark matter and 
MOND-inspired models is yet unsettled even though the need of massive neu- 
trinos to reconcile MOND with cluster data seem to disfavor this latter ap- 
proach. Besides, all analysis on cosmological still assume the existence of dark 
matter. 

It is important to keep in mind that these are not the only tests of the 
deviation from the standard vlCDM that can be performed. Let us mention 

— Test of the weak equivalence principle. They can be performed on a large 
band of redshifts, up to BBN time, by constraining the time variation of 
fundamental constants [6]. 

— Test of the distance duality relation. In standard cosmology the angular 
and luminosity distances are related by Dl = (1 + z) 2 Da(z). This equa- 
tion holds in any metric theory of gravity if the number of photons is 
conserved. By testing it [128j , one can check the validity of Maxwell theory 
and constrain models such a photon-axion oscillation. 

— Test of the Copernican principle. All the equations and solutions we have 
used, assumed the existence of a homogeneous and isotropic background 
spacetime. It is only an assumption based on the Copernican principle and 
recently many proposals to test it appear in the literature [1291113011131] . 

— Propagation of gravity waves. In bi- metric theories of gravity, gravity waves 
and photons may not necessary follow the geodesies of the same metric so 
that there can exist a time delay between them [132 . Confronting their 
arrival times (as well as those of neutrinos, if massive) allows to set con- 
straints on bimetric theories of gravity. If gravity propagates slower than 
light then some tight constraint can arise from the energy loss of cosmic rays 
by gravitational Cerenkov radiation [133] leading to cq^/c— 1 < 2 x 10~ 19 . 

These tests will enable to check the robustness of the hypothesis on which 
our cosmological model rests. It will either confirm the need for the existence 
of dark matter and dark energy (thus extending drastically the domain of 



26 



validity of GR) or offer new theoretical constructions to explain the late time 
acceleration of the cosmic expansion. 
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